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Abstract 

We perform a detailed analysis of the Two-Higgs Doublet Model (2HDM) potential. At the tree- 
level, the potential may accommodate more than one minima, one of them being the electroweak 
(EW) minimum where the universe lives. The parameter space allowed after the data from the Large 
Hadron Collider (LHC) came in almost excludes those cases where the EW vacuum is shallower than 
the second minimum. We extend the analysis by including terms in the 2HDM potential that break 
the Z 2 symmetry of the potential by dimension-4 operators and show that the conclusions remain 
unchanged. Furthermore, a one-loop analysis of the potential is performed for both cases, namely, 
where the Z 2 symmetry of the potential is broken by dimension-2 or dimension-4 operators. For 
quantitative analysis, we show our results for the Type-H 2HDM, qualitative results remaining the 
same for other 2HDMs. We hnd that the nature of the vacua from the tree-level analysis does not 
change; the EW vacuum still remains deeper. 

PACS no.: 12.60.Fr, 14.80.Ec 


1 Introduction 

Two-Higgs Doublet Models (2HDM) [1, 2] are one of the most popular extensions of the Standard Model 
(SM), even without invoking supersymmetry, for which more than one scalar doublet is a necessary 
condition. In a 2HDM with two doublets and 4>2, there are hve physical scalars: the two CP-even 
neutrals h and H, the CP-odd neutral A and two charged scalars While a generic 2HDM may 
contain flavor-changing neutral current (FCNC) interactions (see, e.g., Ref. [4, 5, 6] for such type of 
2HDMs), one usually invokes some discrete symmetry to banish such FCNC at the tree-level, based on 
the Glashow-Weinberg-Paschos (GWP) [7] theorem that there is no tree-level FCNC if all right-handed 
fermions of a given electric charge couple to only one of the doublets. It turns out that there are four 
2HDMs that satisfy the GWP criterion when a discrete Z 2 symmetry is applied on the Lagrangian. 
They are: 

1. Type I, for which all fermions couple with 4>2 and none with <hi; 

2. Type II, for which up-type quarks couple to <I> 2 , down-type quarks and charged leptons couple to 
^1; 

3. Type Y (sometimes called Type III or Flipped), for which up-type quarks and charged leptons 
couple to <I >2 and down-type quarks couple to <I>i; 

4. Type X (sometimes called Type IV or Lepton-specific), for which all charged leptons couple to 
<I>i and all quarks couple to <I> 2 . 

Among them. Type II 2HDM has been most widely investigated [8, 9, 10, 11, 12] because the scalar 
sector of minimal supersymmetry is a Type II 2HDM. 

The tightest constraint on any 2HDM comes from the fact that the observed scalar resonance at the 
Large Hadron Collider (LHC) can be identified with the SM Higgs boson with muiggs = 125.09 ± 0.24 

* indrani300888@gmail.com 
Anirban.kundu.cu@gmail.com 

^ The neutral scalars do not have any definite CP property if the scalar potential violates CP. For 2HDM with 
spontaneous CP violation, see, e.g., Ref. [3]. 
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GeV [13]. Thus, the Yukawa and gauge couplings of the 2HDM must be so aligned as to make the 
lighter CP-even mass eigenstate h to almost coincide with the SM Higgs boson. This is known as the 
alignment limit, and while the allowed parameter space may vary from one 2HDM to the other, the 
qualitative results are quite similar [14, 15, 16]. There are other constraints, like the oblique parameters 
or the lower limit on the mass coming from 6 — 7-57 decay rate [17, 18], but such constraints are 
in general not equally valid for all 2HDMs The theoretical constraints include, just like any other 
extensions of the SM, the vacuum stability, validity of the perturbative nature of the couplings, and 
constraints coming from the requirement of unitarity of scattering amplitudes [19]. 

In this paper, we will focus upon the scalar potential of the 2HDMs. This is much more complicated, 
even at the tree-level, compared to the SM scalar potential, if both the CP-even neutral scalars are 
allowed to have nonzero vacuum expectation values (VEV). As has been shown in Refs. [20, 21, 22, 23, 
24], the 2HDM scalar potential, even with the softly broken Z 2 -symmetry, may allow more than one 
minima. A similar analysis was very recently done for Z 2 -breaking models [25]. In Refs. [20, 22, 23, 24], 
the authors considered the scalar potential of a Z 2 -conserving 2HDM at the tree-level, and showed that 
it can allow multiple normal non-equivalent stationary points (at most two of them can be minima). 
However, a charge-violating or CP-violating minimum cannot coexist with a normal minimum, i.e., 
where the CP-even neutral fields get the VEV. This can be put in a more succinct way; minima with 
different natures cannot co-exist in 2HDM [20, 24]. Also, the data from LHC all but rules out those 
points where the second minimum is deeper than the electroweak (EW) minimum, i.e., the minimum 
where the universe lives. 

In this paper, we consider both Z 2 -conserving and Z 2 -violating 2HDM scalar potentials for our anal¬ 
ysis. By Z 2 -breaking 2HDM, we mean those with dimension-4 operators violating Z 2 , like (<k|<ki)(<I>| 4 > 2 ) 
or (<k 2 ‘h 2 )(‘k|<l> 2 ). There can be soft Z 2 -violating terms like <k|‘k 2 or <k 2 d>i in the potential; if there 
are no dimension-4 Z 2 -breaking terms, we call those models Z 2 -conserving, although, strictly speaking, 
they are not. We show that the conclusions drawn about the nature of the minima of the potential in 
the context of Z 2 -conserving 2HDMs remain valid for Z 2 -violating 2HDMs too. 

To check the robustness of the tree-level results, we further perform a one-loop analysis of the 2HDM 
potential, and choose only those models that show a double minima. A nice review of the one-loop 
corrections in the context of Z 2 -symmetric 2HDM and scale invariant 2HDM can be found in Refs. 
[26, 27, 28]. It is a common knowledge that the one-loop corrections [29] can be significant only in 
the flat direction of the potential. Eor the SM, this is easy to obtain [30]; so is for the Inert Doublet 
models ^ [32, 33] where one of the VEVs is zero. For a generic 2HDM, this is a cumbersome task, 
but can be done, in principle, following the prescriptions of Gildener and Weinberg [34], and the ray 
in the potential space along which the tree-level potential is zero can be found. However, If U 2 ui, 
where vi and V 2 are the VEVs of the CP-even neutral components of $1 and 4*2 respectively (so that 
tan/I = V 2 /V 1 1 ), the potential along the direction is almost flat, so it is instructive to show the 
variations of the potential perpendicular to this direction, i.e., along 4 > 2 . Another important point is 
the setting of the regularization scale ja for the one-loop corrections. Variation of fi is equivalent to the 
variations of the tree-level quartic couplings Xi, as can be seen from a renormalization group argument. 
As the nature of the potential can best be described by these quartic couplings, we would like to forward 
a prescription of choosing /x for the 2HDM: choose it so that the position of the EW minimum remains 
unaltered. As we will see, this keeps the position of the second minimum too almost unaltered. Of 
course, the depths of the potential at the two minima will change. With such a prescription for choosing 
/X, the conclusions about the stability of the EW minimum that were drawn from a tree-level analysis 
remain unaltered. 

^ For example, the & —>■ 57 constraint, nijj± > 316 GeV, can be evaded in Type I and Lepton-specific 2HDM. 

^ For the Inert Doublet model (which is important from the cosmological implication of providing a cold dark matter 
candidate), one of the VEVs is zero, and so it is easier to treat analytically. We, however, will not go into any detailed 
study of such cosmological implications of the 2HDMs in this paper. Another such implication is the successful first-order 
electroweak phase transition and electroweak baryogenesis, which is discussed in Refs. [31]. 
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The paper is arranged as follows. In section II we briefly review 2HDMs, with softly broken Z 2 
symmetry and without Z 2 symmetry, and list the constraints and the minimization conditions on the 
potential. Section III introduces the one-loop corrected effective potential and modified minimization 
conditions for Type II 2HDM with and without Z 2 symmetry. Our results, for both tree-level and 
one-loop analysis, are shown in Section IV. Section V summarizes the paper. Some calculation details 
and relevant expressions are relegated to the Appendix. 


2 Brief review of 2HDM 


2.1 2HDM with softly broken symmetry 


To start with, let us focus on the most canonical 2HDMs, where the Z 2 symmetry is broken only softly 
by a dimension-2 operator, and all the couplings are real. The notations used here essentially follow 
those in Ref. [1]. Later on, we will introduce both dimension-4 Z 2 -breaking operators as well as complex 
parameters in the scalar potential. 

Let us denote the two doublets, both with hypercharge Y = -|-1, by and $ 2 ; which can be written 
more explicitly as 


$ 


a 



xl \ 

(</>a + irja) / ’ 


a = 1,2. 


( 1 ) 


We further assume the VEVs to be aligned towards the direction of the CP-even neutral field, so that 
{(pi) = vi, {(P 2 ) = V 2 , and we conventionally denote tan/3 = V 2 /V 1 . 

Invoking a Z 2 symmetry $1 —)■ 4>i, <I >2 —?• —4*2, so that there is no tree-level flavor-changing neutral 
current (FCNC), one may write 


V (4*1, 4 * 2 ) = mf^4>|4>i-|-77122*^24*2 — m^2 + 4>24 >i^ 


+ ^Ai + ^A2 ($^4*2)% A3 ( 4 *|$i) (^<^ 1^2 

4*1^2)% 


+A4 ($1^2) ($2^1) + ^As 


( 2 ) 


Here mfi, 777 , 22 ) "^ 12 > Ai, A 2 , A 3 , A 4 and A 5 are all real and 777^2 softly breaks the Z 2 symmetry. As 
mentioned before, to differentiate from the potential that breaks the Z 2 symmetry with dimension- 
4 operators, such models will henceforth be called the Z 2 -conserving or Z 2 -symmetric models, even 
though it is broken softly by 777 ^ 2 - 

The two CP-even neutral states (pi and (p 2 are in general not mass eigenstates. The corresponding 
mass matrix can be diagonalized through a rotation by an angle a, and the mass eigenstates are 


h = (p 2 cos a — (pi sin a , H = (p 2 sin a + (pi cos a , (3) 

where h(ff) is the lighter (heavier) eigenstate. 

Note that if \a — j3\ is an odd (even) multiple of 7 r/ 2 , h{H) becomes identical with the SM Higgs 
boson, with a VEV of u = ss 246 GeV. For example, the hVV* (HVV*) coupling is just the 

SM coupling times sin(Q; —/3) (cos(a —/3)), where V is any weak gauge boson, W or Z. The limit where 
h behaves as the SM Higgs boson is known as the alignment limit. The LHC data strongly favours the 
alignment limit and this sets a nontrivial constraint on the parameter space. The allowed parameter 
space, of course, depends on what type of 2HDM is chosen. We refer the reader to Ref. [14, 15] for a 
study of the alignment limit in Type I and Type H 2HDMs, and to Ref. [ 6 ] for a typical example of 
constraints coming from 2HDMs with tree-level FCNC. 
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The most generic Yukawa interactions for these four models can be written as [1], 


Cy = - 


Y, \Y^QLdR^, + + Y^lLlR^j + h.c. 


where = ir 2 <h*, Qr, Lr, cIr, ur and Ir are generic doublet quarks, doublet leptons, singlet down-type 
and singlet up-type quarks, and singlet charged leptons respectively. Y^, Y? are 3x3 complex 
matrices, containing Yukawa couplings for the down, up, and leptonic sectors respectively. In our 
analysis we will consider only top, bottom, and r Yukawa couplings to be nonzero. 

The masses of the charged Higgs, H^, and the pseudoscalar. A, can be written as 


mH± = 


m\ = 


^ (A4 + As) , 
V 1 V 2 2 

— Asu^ . 

V 1 V 2 


The couplings can be expressed in terms of masses of the physical states and the mixing angles a 
and P as [ 20 ] 

\ ^ 2 I 22 


Ai = + 

\ ^[ 22,22 


A 3 — 


/ 2 _ 2 \ _ ^12 

v"^ v‘^S 2I3 ^ ^ v'^SpCp 


^4 = -2 - 2m^±) + 




As — 


^12 _ ^ 
v'^SjiCp 


where mh{mH) is the mass of h{H), and cq and sg are generic shorthand notations for cos0 and sin0 
respectively. 

The requirement that the scalar potential always remains bounded from below leads to the following 
stability conditions for 2HDMs with Z 2 symmetry [1], 

Ai,A 2 > 0 , A 3 > —\/AiA 2 , A 3 -|- A 4 — |As| > —\/AiA 2 . (7) 

Note that A 3 , A 4 , and As can in principle be negative; however, if the Z 2 symmetry is exact, this may 
lead to tachyonic masses for and A. 


2.2 2HDM with hard Z 2 breaking terms 

The most general renormalizable scalar potential of 2HDM can be written as [1] 

Y(^>i,^>2) = + ml2^l^2 - ml2 (^1^2^ + h.c. 

+ ^Xi -F ^A 2 ($2^2) + As (^2^2) + A4 (<I> 1 ^ 2 ) (^>2^1) 

+ ^As (^1^2)^ + Ag (^1^2) + At (^^^2) (^1^2) + h.c. , ( 8 ) 

where the model parameters rnh, m| 2 , Ai, A 2 , A 3 , A 4 are real and mf 2 , As, Ag and A 7 can in principle 
be complex, and “h.c.” stands for hermitian conjugation. 

We will follow the same convention as the previous subsection and consider two different cases, 
namely, mf 2 , A 5 , Ag and A 7 are (i) real, and (ii) complex. 
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2.3 2HDMs without Z 2 symmetry with all real parameters 


For this case, the masses of and A can be written in an analogous way of Eq. (5), using the 
minimization conditions of the potential: 

2 1 2 2 ^ \ \ \ 2 \ 2 ^ 2.2 

m^± = - mf 2 V - 7 ; (A 4 + As) - —Xev - —Mv , 

V 1 V 2 2 2 v 2 2 vi 

2 1 2 2 \ 2 \ 2 '^2 . 2 7n^ 

= - m^ 2 ^ - A 5 V - —Xev - —Xjv . (9) 

V 1 V 2 2 v 2 2 vi 

The requirement that the scalar potential always remains bounded from below leads to the same set of 
equations as in (7), with an extra condition: 


Ai,A2>0, As > —\/ A1A2 , As + A4 — jAs] > —\/ A1A2 , 

2|A6 + A 7 I < - (Ai + A 2 ) + As + A 4 + As . (10) 

The mixing angle between the CP-even neutral states is given by a = ^ arctan(yl/i3), where 


A 

B 


3 3 

^-12 “ (A 3 + A 4 -I- As) V 1 V 2 — 2 ^ 6 '^! — 2^7V2 

-^AiUi -h ^A2 u| -h {vi/v2 - V2/V1) + ^ (A 7 - Ae) viV2 + ^A 7 uf - ^Aeu? . (11) 


2.3.1 Minima of the potential 


The scalar potential of the 2HDM shows a much more complicated structure than that of the SM. The 
potential can have multiple non-equivalent normal stationary points (only two of them can be minima) 
[20, 22, 23, 24], depending on the parameters. These minima can be all normal (where only (pi and (p 2 
get nonzero VEV), charge-breaking (where at least one of the charged fields xt S^ts a nonzero VEV) 
or CP violating (where the two VEVs have a nontrivial phase between them). It has been shown in 
Refs. [20, 21, 22, 24], that for the canonical Z 2 -symmetric 2HDM, (i) existence of a normal minimum 
rules out a charge breaking or CP violating minimum, those extrema can only be saddle points at 
best; (ii) there can be more than one normal minima, only one of which corresponds to the Standard 
model (SM), which we call the EW minimum. If the EW minimum is the global one, we are in a stable 
situation; if it is not, the universe may tunnel down to the deeper minimum if the tunneling time is less 
than or of the order of the lifetime of the universe. It has also been shown that the LHC data effectively 
rules out the parameter space where the EW minimum is shallower. We would like to extend this result 
to Z 2 -breaking 2HDM with both real and complex parameters, and also investigate the nature of the 
potential when one-loop corrections are taken into account. 

In a normal minimum, we may write 


{^i)n 


1 



, {^ 2 )n 


1 



( 12 ) 


In certain situations, several non-equivalent normal extrema are allowed by the minimization conditions 
of the potentials and at most two of them can be local minima. Thus, there can be a second minimum 

with VEVs {v[,V 2 ), with v' = \Jv'\ -|- v '2 7 ^ 246 GeV, where exactly same symmetries are broken. 
Obviously, even with the same parameters of the potential, the masses of all SM particles are going to 
change from the EW vacuum configuration. 

At the EW vacuum, the potential, as follows from Eq.( 8 ), can be written as 


Vo 


—ni\iv\ + —7n^2'^2 ~ 7 ni 2 ViV 2 + g Aiu^ -|- -X 2 V 2 + - (A 3 -|- A 4 -|- A 5 ) 
+ ]^XqvIv 2 + ]^X^vlvl , 


( 13 ) 
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and at the other minimum, we get the corresponding Vq by replacing {vi,V 2 ) in Eq. (13) by (^'^,^ 2 ). 
If Vo < (>)^o^ (second) vacuum is stable and is the global minimum. 

The minimization conditions of the potential are easy to obtain: 


fl{vi,V2) 

f 2 {vi,V 2 ) 


m\iVi — rni2V2 + 2 + -^4 + ^5) V1V2 + 

^22'^2 ~ '^'12^1 + 2'^2'i^2 + 2 + A4 + A5) U2U4 + 


3 1 

-X 6 vfv 2 + -X 7 V 2 

1 3 

-Aguf + -Xjvlvi 


0 , 

0. (14) 


with 


3 1 

mil 4“ 2 '^^^! 2 + ^4 + A 5 ) + 3 A 6 U 1 U 2 

3 1 

^22 4" 2^^^2 4" 2 (^3 4- A 4 + A 5 ) vf + 3 A 7 U 1 U 2 


> 0 , 
> 0 . 


(15) 


2.3.2 Charge-breaking minimum 


If the symmetry breaks in such a way that one of the charged helds, say X 2 > ^ non-zero VEV 7 

too, such that 


the minimization conditions of the potential are 




2)CB 


1 

71 



(16) 


and 


mjivi 

m 22 V 2 
17122 4 “ 


- mj2V2 + ^AiU? ^AsUi {V2 + 7 ^) + ^ (^4 + A5) V1V2 + ^A6UiU2 + ^At (7^ + V2)v2 = 0 , 

— mi2Vi + -X2 (u| -|- 7^) U2 + - (A3 -|- A4 -|- A5) V2vf + 2'^6'yi 4- -A7 (7^ -|- 3 u|) ui = 0 , 

^A2 (7^ + U2) + ^Asuf-h A7U1U2 = 0 , ( 17 ) 


3 1 1 

mil 4“ 2^3 ^i 4“ 2^^ (4'^ 4- ul) + - (A4 -|- A5) + 3A6UiU2 > 0 , 

^22 4" 2-^2 (3u| -|- 7^) -|- - (A3 -|- A4 -|- A5) vf + 3X7V1V2 > 0 , 
^22 4- ^A2 (37^ -I-u|) + ^A3 Ui -I-A7U1U2 > 0. 


(18) 


2.3.3 CP violating minimum 

At the CP violating extremum, one of the CP-odd neutral fields acquire a nonzero VEV S. In other 
words, the VEVs develop a relative phase. Thus, 

((</'! + ^'l) + i ivi + <5))) ’ iih + V 2 ) + im) 

We have similar minimization conditions for the potential: 

m\iVi — m‘l2‘V2 4 - —Ai (vf -|- 2 4 - A4 -|- A5) viv\ -\- -Ag (3^4 -|- < 5 ^) V2 + —X7V2 = 0 , 

m22V2 ~ m‘l2‘Vi + 2^2'1’2 4 “ 2 (^3 4 - A4) V2 (vf -|- 5 ^) -|- —A5 (^4 — ( 5 ^) V2 + —Ae (vf + 5 ^) U4 

3 . 2 . 

4--A7U1U2 — 0, 

mfi6 + ^Xi ( 5 ^ + vf) 5 +]^ (As -h A4) vfS - ^Asufd -h XqViV25 = 0 , ( 20 ) 
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and 


mil + 2 ^^ (^'^1 + + 2 + ^4 + As) V 2 + 3A6Wi'V2 > 0 , 

3 1 1 

"^22 + 2 2 ”^^ ~ SX7ViV2 > 0 , 

mfi + —Ai ( 3 i ;4 + (5^) + — (A 3 + A 4 — A 5 ) V 2 + A 6 Wif 2 > 0. ( 21 ) 

2.4 2HDM without Z 2 symmetry with complex parameters 

Next we consider the 2HDMs with 171 ^ 2 , A 5 , Ae and A 7 complex, and use 

mi 2 = ml 2 + iiTT-ii > -^5 = A 51 + iA52 , Ae = Aei + 1 X 52 , Xj = A 71 + iXj 2 ■ (22) 

Most of the relevant expressions are identical with the real parameter case discussed before, with 
the (mf^, A 51 , Aei, A 71 ) replacing (mf 2 , A 5 , Ae, A 7 ) respectively. While this applies to the expressions 
of physical mass eigenstates, there is an extra condition so that the charged Goldstone boson is still 
massless: 

— 77732 + 2 A527’17’2 + 2 '^ 62 'ri + ~X72V2 = 0 . (23) 

This shows that only three out of the four phases are independent. 

The conditions for the stability of the potential are 

Ai, A 2 > 0 , A 3 > — y^AiA 2 , A 3 + A 4 — ij AI 3 + A 52 > —\/XiX 2 , 

2|A6i + A 71 I < - (Ai + A 2 ) + A 3 + A 4 + A 51 , 2|A62 + A 72 I < - (Ai + A 2 ) + A 3 + A 4 — A 51 . (24) 

Similarly, the minimization conditions for the normal and charge-breaking minima are obtained 
from the corresponding expressions for the real parameter case by the substitutions mentioned before. 
For the charge-breaking minima, the imaginary parts of the potential parameters play a role. The 
minimization conditions are 

mjivi - m\ 2 V 2 + ^Ai [vl + 6 ^) ^ (A 3 + A 4 -k A 51 ) -k ^A52u|(5 

T^Aei (3 i; 3-k (5^) U 2 -k A62'yi7’25-k ^A7 iu| = 0 , 

m 22 'y 2 - mf^vi - ml 26 + ]^X 2 vl + ^ (A 3 -k A 4 ) V 2 {vj + 5^) -k ^Asi [vj - V 2 

1 1 3 3 

+A52'yi'y2<^ + 2"^61 (^1 ^ = 0 ) 

mfiS - mliv 2 + ^Ai {S^ + t’l) ^ (A 3 + A 4 ) v^d - ^Asiu^h -k ^A 52 Uiu| 

+A61'yi7'25 -k 2'^62 (3h^ -k U 3 ) V 2 + —X 72 V 2 = 0 , (25) 

and 

mfi + -Ai ( 3 U 3 -k (5^) -k - (A 3 -k A 4 -k A 51 ) -k 3 A 61 U 1 U 2 -k A627’2<^ > 0 , 

3 1 1 

"^22 + 2 ^ 2 vI -k - (A 3 -k A 4 ) {vf -k (5^) -k -A 51 {vf - (5^) -k A 52 U 15 -k 3 A 71 U 1 U 2 -k 3 X 72 V 2 S > 0, 

mfi + —Ai (Svf + (5^) -k — (A 3 -k A 4 — A 51 ) V 2 -k A 61 U 1 U 2 -k 3 Xq 2 ‘V 2 ^ > 0. (26) 

The tree-level analysis will be performed using these expressions and the corresponding solutions of 
the simultaneous minimization equations. 
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3 One-loop effective potential 

The form of the one-loop correction to the scalar potential is quite standard, and is given by 


mi mi{(l)ci,4>c2), Vi = ^ Nimf ^In ^ - C* j , (27) 

where the sum runs over all bosonic (the physical scalars, the unphysical Goldstone bosons, W and 
Z) and fermionic (for our case, only t, b and r) degrees of freedom. Here n is the regularization scale. 
The masses are field-dependent quantities, being functions of the background fields (f)ci and (j)c 2 - They 
can be thought as the positions of the minima in the field space. In the limit 4>ci = vi, 4>c2 = V 2 , the 
field-dependent masses are equal to the physical masses. The constants Ni and Ci are given by 

Nh = Nh = Na = Nqo = 1 , Nh± = Nq± = 2 , 

Nw = Q , Nz = 3, Nt = Nb =—12 , Nr =—4:, 

Ch = CH = CA = Ch± = Cgo = Cg± = a = a = a = ^, Cw = Cz = \. (28) 

The full potential can be written as 

V = V{<l>u^2) + Vi (29) 

where H(<I>i,<I> 2 ) is the tree-level potential as given in Eq. (2) or (8). 

At this point, let us again note that the masses in Eq. (27) are functions of background fields (/>ci 
and (f>c 2 - The positions of the minima are functions of not only (/>ci and (j)c 2 but also the regularization 
scale fi. To have any idea of the nature of the potential after one-loop correction, one has to fix // 
by some prescription. The renormalization group improved one-loop potential clearly shows that the 
variation in ^ is equivalent to the redefinition of the coupling parameters of the theory. One popular 
way is to fix it in such a way that the one-loop corrections are minimum, hoping that this will minimize 
the higher-order corrections too. To be physically more transparent, we try a different approach: we 
fix n in such a way that the position of the EW minimum remains (almost) unchanged with respect to 
the tree-level position. This will keep the one-loop corrected field-dependent masses to be at the same 
values of the tree-level masses at the EW vacuum. Only the depth of the potential changes by the 
one-loop correction. Thus, our prescription is to tune fj, in such a way that (pd = vi and (j)c 2 = V 2 (so 
that, for example, the Goldstone bosons are still massless ^). As we will show later, if one tunes /x so 
that the position of the EW minimum is unchanged, it will keep the position of the second minimum 
almost unchanged too. 

As a concrete example, let us now discuss the Type II 2HDM; as the Yukawa couplings enter 
the picture, one needs to specify the type of 2HDM under consideration. We have checked that the 
qualitative features remain unchanged in all other 2HDMs. 


3.1 Type-II 2HDM without Z 2 : One-loop corrected potential 

To keep the discussion as much general as possible, let us focus on the one-loop correction to the 2HDM 
(without Z 2 symmetry) tree level potential. We can obtain the same for Z 2 symmetric 2HDM by putting 
Ae and A 7 to be equal to zero, and by making A 5 and 771^2 real. Using the definitions of fi{vi,V 2 ) from 
Eq. (14), we can write the modified minimization condition for one-loop corrected potential V as 


fi{vi,V2) 


dVi 

dcpi 


= 0 , 


f2{vi,V2) 


dVi 

d(p 2 


= 0 , 

</>c2='f2 


(30) 


■^The treatment of the Goldstone bosons in one-loop corrected potentials is tricky, and a consistent treatment needs 
resummation of the Goldstone contributions in the effective potential [35]. 



where, 

dVi 
d(f)i 


cl=Vl 


647r^ 


Am\F{h, 1) ^In ^ ~ Am%F{H, 1) ^In - 1^ 


m\ 


+4m^ 1^1 - In l^'ViAsi - - 

+4m|-± ^1 - In-^^ ^'ViA 45 


2 R 

12 


m 


1 - -4 


V 2 


-egjm^vi 


1 , m 

-In — 




- Sg'^mzvi 


+8Y^m‘^vi ( 1 — In ^ 




4v2 


2 v 2 


2 R 

12 


1 - - 

-L O 


+ -^61^^ Aeit^f 


2112 


V2 


X7lV2v‘^ A7i1;2\ 

) 

X 7 lV 2 v'^ 

- 2 - 1 " A71^^2 


2 i; 


m7 




1 

mf 


(31) 


and 


dVi 

d4>2 


4>c2='V2 


6471^ 


mt 


+^m\ 


1 — In 


m, 




g- 


V 2 X 51 - 


AmiF{h, 2 ) { In ^ - 1 + AmijF{H, 2 ) In ^ - 1 


m 


g^ 


m 


2 R 

12 


2 i;i 


1 - 


Vn 


Xii'v'^ Xviv- 

+ 4 — + 


+4m^± ^1 - In ^'V2A45 - 


^2R 

^12 

Vi 


i-d 

vl 


4vi 2i;i 
+ ^71^'2 


Xeiviv^ AeiwA 

2 J 


2 v\ 


Vl 


4vl 

Aei'^i^’^ , , 
-2-^ Aeifi 


2 i;o 


- 6 gjmwV 2 


i-ln^ 

3 g^ 


- 3g^mzV2 ( ^ - In + 24y/mj 112 Tl - In ^ 

6 g j \ g 


(32) 


with A 45 = A 4 + A 51 and g = \/g\ F g\, gi and g 2 being the U{1)y and SU{2)l gauge couplings. Yukawa 
couplings for t, b, and r are denoted by Yt, Y^ and Y^ respectively. The F-functions have been defined 
in Appendix A. Note that these expressions are valid only if we tune g to keep the position of the EW 
minimum of the one-loop corrected effective potential V, defined in Eq. (29), unchanged with respect 
to its position in the tree level potential. One also needs to check the second derivatives to ensure that 
the extremum is a local minimum: 


d(j)j 


> 0 , 

4>cl=Vl 


d(t>l 


> 0. 


<t>c2=V2 


(33) 


The expressions for second derivatives are also given in Appendix A. Note the absence of the Goldstone 
bosons in Eqs. (30,31,32) because of our choice of g which keeps them massless. 


4 Analysis and Results 

4.1 2HDM at tree-level 

Erom a random scan over the parameter space spanning over 7 x 10® different choices of model param¬ 
eters, we generate a number of models for which the following conditions are satisfied. Our analysis 
includes the canonical 2HDM with Z 2 symmetry, and 2HDM without Z 2 with both real and complex 
parameters. 

• The potential has to be stable at all scales before it either blows up (due to one or more couplings 
hitting the Landau pole) or becomes unbounded from below. 

• The dimensionless couplings must remain perturbative over the entire range of validity of the 
theory, except maybe at the very end where they approach the Landau pole. 
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• There should be one minimum of the scalar potential for which v = \Jv\^v\ = 246 GeV. This 
we will call the EW vacuum. This, in fact, acts as the tightest constraint and rules out the largest 
chunk of randomly generated models. We focus only on those models that allow another local 
minimum apart from the EW vacuum. Both tree-level minimization conditions, for vi and V 2 , 
should be satisfied in both the vacua. 

• Other constraints like that on the charged Higgs boson mass coming from 6 — 7-57 satisfied. 
We use mfj± > 316 GeV. While not all constraints are valid for all 2HDMs, we focus only on the 
Type II 2HDM. 

• The 125 GeV resonance found at the Large Hadron Gollider must have properties close to that 
of the SM Higgs boson. In other words, the alignment limit should be maintained. We have kept 
\a — /3| to be between 0 . 97 r /2 and l.l 7 r/ 2 , noting that this is a rather conservative limit. 

The ranges for our scan is as follows: 

0.0 < Ai,A 2 < 1.0, -1.0 < A 3 ,A 4 ,A 5 < 1.0, m? 2 < 4 xl 0 ® GeV^ 1.0 < tan/?< 50.0 , (34) 

where the parameters are taken to be at the EW scale. 

Apart from this, for the Z 2 -violating cases, we have taken 

— 1.0 < Ae, Ay < 1.0 (real couplings), — 1.0 < A 51 , A 52 , Aei, A 62 , A 71 , A 72 < 1.0 (complex couplings). 

(35) 

Also, for this case the scan is made on 771^2 = ^12 + as the phase in mf 2 is fixed by Eq. (23). 

We find a few common characteristics for these models that allow two minima. They are: 

(i) The EW vacuum is always deeper than the other vacuum. This happens mostly because of the 
imposition of the experimental constraints. Thus, even with the introduction of Z 2 -breaking parameters, 
there is no chance of tunneling to the other minimum, at least at the tree-level. This reinforces the 
conclusions obtained by the authors of Ref. [20]. 

(ii) If both the vacua are normal, there is no other minimum that breaks charge conservation or CP. 
This, again, is in tune of what the authors of Refs. [22, 24] found. 

4.2 2HDM at one-loop level 

We would now like to perform the same analysis on the one-loop corrected potential on those models 
that satisfy the initial constraints and show the presence of a double minima at the tree-level. The 
regularization scale fi is so chosen as to make (/>ci = vi and (/)c 2 = V 2 . This is, of course, a rather 
restrictive choice, but keeps all the masses as well as tan/3 unchanged at the EW minimum even after 
the one-loop corrections are implemented. 

In a generic 2HDM where both CP-even neutral fields can get nonzero VEV, this is a tricky job, 
and mathematically much more complicated than the cases of only SM, or SM extended by a gauge 
singlet scalar, or the inert doublet models where one VEV is always zero. The complication is further 
enhanced by the fact that at the tree-level, there are two minima of the potential. 

What we do is the following. We fix (fid = vi and tune ^ in such a way that (pc 2 coincides with 
V 2 . This keeps the position of the EW minimum invariant but changes the depth. We could have done 
this the other way round too, namely, keeping (j)c 2 = V 2 and adjusting ^ to make (/)ci coincide with 
the tree level value; however, we prefer the first approach as the contribution of V 2 is larger in v for 
tan/3 > 1. We then use the same ^ but fix the classical minimum for (f>i at ® (p'd = v'l- It so happens 
that (ji'd, in all cases, lies close to the coincidence is not as exact as the EW vacuum, but this being 
the shallower secondary minimum, we will not be so much bothered about its exact position. Thus, we 

®We use primes to denote the corresponding qnantities in the second minimum. 
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compare F(/U, (/>ci = vi, (j)c 2 = "^ 2 ) with = v'^, (p '^2 ~ '^ 2 )- We will also show the effect of varying 

the scale /x on the classical field values. 

What we observe is that the conclusions drawn from a tree-level analysis is more or less unchanged; 
the deeper minimum remains deeper. This is not entirely unexpected, as the one-loop corrections are 
only a small effect that cannot overcome the difference in depths of these two minima. For tan /3 ^ 1, 
4 >i direction is almost flat, so we show our results in the constant-(/)i direction, varying the potential V 
with ( 1)2 ■ 

Our results are discussed in more detail later for some typical benchmark points, but for both the 
cases that we consider {Z 2 conserving and breaking), we never found the second minimum becoming 
deeper than the EW one after the one-loop corrections. Thus, if the EW vacuum is the deeper one 
at tree-level, it remains so after the radiative corrections; there is no chance of developing a deeper 
vacuum and tunnelling into it. In fact, if /x is tuned in such a way that the EW vacuum is not shifted 
from its tree-level position, it always gets deeper by the one-loop corrections] Vi at the EW vacuum is 
always negative. 

4.3 One-loop corrected 2HDM with Z 2 symmetry 

In Table 1, we show three benchmark points for the Z 2 -conserving Type-II 2HDM, characterized by 
small, medium, and large values of tan/3(= V 2 /V 1 ) respectively. All these models show the existence of 
a second and shallower minimum compared to the EW one. These three models, namely, .^2^1, .^2^2, 
and Z 2 C 3 are valid up to 3.8 x 10^ GeV, 2.8 x 10^^ GeV, and the Planck scale respectively; for the 
first two models, at least one of the couplings become nonperturbative at the validity scale and the 
model soon hits the Landau pole thereafter. Obviously, the comparatively low validity range for the 
first benchmark can be ascribed to the relatively large quadratic couplings to start with at the EW 
scale. The evolutions are checked with one-loop renormalization group equations for the Type II 2HDM 
[1, 36]. 



Benchmark 

Parameter 

^2^1 

Z 2 C 2 

Z 2 C 3 

Ai 

0.413 

0.642 

0.068 

A2 

0.842 

0.328 

0.260 

A3 

-0.265 

0.065 

0.132 

A4 

-0.720 

-0.365 

-0.489 

As 

-0.929 

-0.786 

-0.498 

mfi (GeV^) 

2.3 X 10® 

6.2 X 10® 

2.61 X 10® 

ml2 (GeV2) 

2.8 X 10® 

5.2 X 10^ 

-7.77 X 10^ 

mf 2 (GeV^) 

8.10 X 10® 

9.11 X 10^ 

4.84 X 10^ 

vi (GeV) 

83.76 

37.65 

5.05 

V 2 (GeV) 

231.30 

243.11 

245.95 

v[ (GeV) 

442.46 

452.48 

399.18 

v '2 (GeV) 

872.82 

953.66 

775.94 


Table 1; Benchmark points for Z 2 -conserving Type-II 2HDM. 

The positions of the one-loop corrected minima and the corresponding regularization scales are 
shown in Table 2. The benchmarks are chosen scanning the range of tan (3 as well as the other param¬ 
eters. 

The potential profiles are shown in Eigure I. These are drawn as a section of the actual three- 
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tan 13 

(pel 

<Pc2 


^cl 

(t>'c2 

Z 2 C 1 

2.76 

83.76 

231.30 

976.0 

442.46 

846.7 

Z 2 C 2 

6.46 

37.65 

243.11 

492.5 

452.48 

932.6 

Z 2 C 3 

48.7 

5.05 

245.95 

335.0 

953.66 

774.6 


Table 2: The one-loop corrected minima (all quantities except tan/3 are in GeV) for the three bench¬ 
marks. Note the tiny shift of (j )'^2 from t> 2 - 


dimensional plots, for fixed values of 4>ci- That is why the second minimum is not apparent; it occurs 
at a different value of (pci- Also, the tree-level potential does not go to zero as </>c 2 —>• 0, unless (pci 
is tiny, as in Z 2 C 3 . Note that the one-loop corrected potential is always deeper than the tree-level 
potential at the EW minimum; Table 2 shows that the second minimum (p '^2 almost coincides with 





Figure 1: The plots of the tree-level and one-loop corrected potential, with the section taken at a fixed 
value of (pel as indicated on the plots and in Table 2. The upper panel plots are for benchmarks Z 2 CI 
(L) and Z 2 C 2 (R), while the lower panel plot is for Z 2 C 3 . In every plot, the upper curve (blue) denotes 
the tree-level potential profile, while the lower one (red) is for the one-loop corrected potential. 

Only if (pci is small, like in Z2C3, the (pi direction can be approximated by a flat direction. The 
flatness is really impressive: for a 10% (1%) change in <pci, the potential changes only by 0.02% 
(2.5 X 10“^%). However, we have not found any case where the one-loop corrections remove the second 
minimum. 

If we keep (pd = vi or and play with /r as a free parameter, (pc 2 and (p '^2 changes. In Fig. 2, we 
show how (pc 2 and (p '^2 change with // for the three benchmarks. 
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Figure 2: Variation of (j)c 2 (L) and (R-) with /r for Z 2 CI (blue), Z 2 C 2 (red), and Z2C2> (golden). 


4.4 One-loop corrected 2HDM without Z 2 symmetry 


The analysis is analogous to what was performed for the Z 2 -symmetric case. The potential profiles 
are shown in Fig. 3. For the three benchmark points .^ 21^15 -^21^2, and Z 2 V 3 , the ^-values are fixed 
at 641.0 GeV, 655.3 GeV, and 2228 GeV respectively. While (pel = 'Wi, (pc 2 = V 2 , and (p'^i = v[ were 
ensured, the (p '^2 values were quite close to that of they are at 506.2 GeV, 708.9 GeV, and 1783 GeV 
respectively. 



100 


200 300 

02lGeV] 


400 


500 


600 




02[GeV] 


Figure 3: The plots of the tree-level and one-loop corrected potential, with the section taken at a fixed 
value of (pci as indicated on the plots. The upper panel plots are for benchmarks .^ 2 1^1 (L) and Z 2 V 2 
(R), while the lower panel plot is for Z 2 V 3 . In every plot, the upper curve (blue) denotes the tree-level 
potential profile, while the lower one (red) is for the one-loop corrected potential. 
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Benchmark 

Parameter 

^2^1 

Z 2 V 2 

Z 2 V 3 

Ai 

0.656 

0.342 

0.497 

A2 

0.188 

0.928 

0.456 

A3 

0.836 

0.998 

0.089 

A4 

0.659 

0.375 

0.598 

A51 

0.764 

-0.956 

-0.533 

A52 

0.163 

0.666 

0.923 

Aei 

0.633 

0.735 

0.680 

A 62 

0.209 

0.619 

-0.929 

A71 

-0.820 

-0.911 

-0.810 

A72 

-0.0016 

0.709 

0.774 

mfi (GeV2) 

7.47 X 10® 

1.06 X 10® 

1.34 X lO’^ 

m^2 (GeV^) 

3.1 X 10^ 

1.78 X 10^ 

3.71 X 10*^ 

mff (GeV^) 

4.92 X 10® 

1.71 X 10® 

7.78 X 10® 

mji (GeV2) 

3.92 X 10® 

2.49 X 10^ 

2.49 X 10^ 

vi (GeV) 

129.13 

43.89 

14.67 

V 2 (GeV) 

209.39 

242.05 

245.56 

v[ (GeV) 

295.20 

257.06 

355.55 

v '2 (GeV) 

631.54 

799.49 

2064.3 


Table 3: Benchmark points for Z 2 -violating Type-II 2HDM. Note that Aj = Aji + i\i 2 for i = 5,6, 7. 
Also, m ^2 = + *^ 12 - 




Figure 4: Variation of (I)c 2 (L) and (R-) ^ Z 2 VI (blue), Z 2 V 2 (red), and Z 2 V 3 (golden). 
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Just like the Z 2 -conserving case, we choose three benchmarks for the case where the potential 
violates the Z 2 and also involves complex parameters. These benchmarks are shown in Table 3. 

The flatness of the potential in the (/>i direction is again most manifest for the benchmark 

with lowest ui. At the same time, the EW minimum always gets deeper with the one-loop corrected 
potential, and there is no qualitative change from the tree-level result. Because of the comparatively 
large values of the couplings, all the three models tend to hit the Landau pole much before the Planck 
scale, namely, at 45.8 TeV, 4.8 TeV, and 216.7 TeV. This is expected because Z^y^ starts with larger 
values of the quartic couplings at the EW scale. In general, for the double-minima case when Z 2 is 
broken, most of the couplings at the EW scale have to be large to start with and thus such class of 
models are not stable beyond a few hundreds of TeV, which may be contrasted with the Z 2 -conserving 
double-minima models. 

In Fig. 4, we show how (^c 2 and i;A ^2 changes with ^ for the three benchmarks. 


5 Summary 

In this paper we have tried to investigate the nature of the potential of Type II 2HDM, breaking the 
Z 2 symmetry either softly or through dimension-4 operators. There are some known results for the 
Z 2 -symmetric 2HDM at the tree level. Our goal was to investigate how far these conslusions are reliable 
if one (i) breaks the Z 2 symmetry at the tree level with operators with real or complex couplings, (ii) 
does a one-loop correction on the potential. 

For the first part, we find that the introduction of Z 2 breaking does not change the conclusions 
qualitatively: the scalar potential can accommodate at most two local minima, both of which have to 
be normal. If there is a normal minimum, there cannot be a charge-breaking or CP violating minimum 
of the potential. The LHC data highly disfavours those models where the second minimum is deeper 
than the EW minimum, possibly making the EW vacuum an unstable or metastable one. 

For one-loop corrections, we use a regularization scale that keeps the position of the EW minimum 
invariant, changing only its depth. When we focus on models with two minima, this prescription keeps 
the position of the second minima almost unchanged too. The one-loop corrections cannot change the 
relative depths of the minima, i.e., the EW minimum still remains deeper after the correction, ruling 
out the possibility of a metastable vacuum. The conclusions are identical for Z 2 symmetric and Z 2 
breaking 2HDM. While the conclusions were drawn for Type II 2HDM, the results are qualitatively 
the same for other 2HDMs too, as the only change comes from the Yukawa couplings that enter the 
one-loop corrections. 

The stability of the 2HDM at higher energy scales is a complex issue because of more fields and 
couplings. One needs the stability conditions to be valid at all scales and the couplings to remain 
perturbative for calculability (or, at least, not blow up). A general tendency that can intuitively be 
deduced from the RG equations is that higher values of quartic couplings at the EW scale tend to 
pull down the range of validity of the theory, which means that some other ultraviolet complete theory 
takes hold beyond that range. However, a large part of the parameter space is still compatible with the 
stability up to the Planck scale for the Z 2 -conserving class of models. For the Z 2 -violating class, the 
existence of two minima generally forces some of the quartic couplings to be large at the EW scale and 
hence such models cease to be valid beyond a few hundreds of TeV at the most; for smaller couplings, 
one gets the single-minimum models. 
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A Expressions for first and second derivatives 

We use the following shorthand notations; 

(A.l) 


F{a, i) = mo 


dma 

dcpi 


G{a, i) = 


i=Vi 


dmc 

d(j)i 


+ rUc 


mo 


d(f): 


Fih,l) = 


4 V ^11^2 ^ ^ 2; 2 2 vj 2 v2 J 

-4_ ( („? - „!) + - A,i) _ ^ 

4 :^/a V viV 2 2 2 ui 2 v 2 J 


/ Yr)^R rfY)‘^R 

Uv,-^ + {vl-vl)^ + 

V V 2 ^ ' 2 v{V 2 


3u 2 (Aei — A71) A71U2 SAeiufA 


1 


2 y/a 


—ml^ + A345Uir;2 + 


2vIv2 ' 4 

SAeiuf 3A71U2 


+ 


4vl 


+ 


4u2 j 


+ 


X (A345'f^2 + 3A61 Ui) , 


Fih,2) = - 


1 , "^12 f 2 2\ , 3 A 67 t’i 3 A 7 iu| Aei^n 

4 ("2 - -4 + ^-^ ^ j 

1 A „,2 \ „,2 7„,2 „,27 , 3 r;iU 2 (Aei - A 71 ) A 7 iui , Aeiuf 

-JTi ^-^ ^ 

f X , "^1^ , / 2 2 N "iff , 3ui(A6i-A7i) 3A7iu| AeiufA 

+ — + («, -I'J ^ +-J- 


1 


2 ^ 


+ A 345 U 1 U 2 + 


3A6iuf 3 A 71 U. 


+ 


X (A345U1 + 3A71U2) , 


mi) = 7 


+ ‘ (Ai 7 - ^ (vf - , 1 ) + - At.) _ ^ ^ 

4^/0 \ uir;2 2 2ui 2^2 / 

/ ™ 2 i? 2 /? 

(A4U4-m + (,2_ 2)^^ 

V V2 ' 2 vfV 2 


3^2 (Aei — A71) A71U2 3 A 6 iufA 


+ 


2 x/a 


—ml^ + A345Uir;2 + 


2vIv2 4 

3 A 6 iuf 3A71U2 


4vl 


4r;2 / 


+ 


X (A 345 U 2 + 3A61Ui) , 


F{H,2) = - 


1 , "^12 (2 2 \ , 3 A 67 t’l 3 A 7 iu| A 6 iuf\ 

4 + il7 ^-^ ^ j 

+J- ( A..; - X-o 4 - ^ (7 - vl) + - An) _ ^ 

4V^ V ^^^2 2 2ui 2 u2 ) 

( X , "^12^ ,72 2A "iff , 3 ui(A 6 i-A 7 i) 3 A 7 iu| AeiufA 

+ — + (^. -I'J ^ +-J- 


+ 


2x/a 


—ml^ + A345Uir;2 + 


3A6iuf 3A7 iu| 


X (A345r;i + 3A7ir;2) , 


(A.2) 
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where 


^67 — Aei + A71, A345 — A3 + A4 + A51, 

^ _ f (Alt’! — X2V2) 3r;ir;2 (Aei — A71) 




4rii 


+ 


mff {vf - vl) 


Av2 


+ ( -^12 + A345'yi^^2 + 


SAeiff SAt^Y 

2 2 ) 


V 1 V 2 


(A.3) 




d(t)\ 


4’cl='Vl 


9Vi 


d4>2 


4 'c 2 ='V 2 


1 

6477^ 

+8 




V 


Vi 


3A61W1 , A7i7;2 
— A 51 -^-h 


m 


,2R 

12 


2i;2 


1 - 1 ) “ ^^lAsi — 

v( 


2V2 
Aeiv^ 


2i;i 

Aei'wf 


A7ii;2f" 

2i;3 

A71'1’2 


V - V X . 4112 2172 2 

, o 2 f mj^V2 1, 3A61W1 , A7ii; 2 A7ii;2i;^ 


+ 


+4 


™2R 

^12 

V2 

2 2 


1 - 1 ) “ ^iA45 — 


2v2 2i;i 

AeiT’i _ A6i17^ 

V2 2v2 


m\ 

In^ 


A7ii;2f‘ 

4i;f 


- 1 


— A7ii;2 + 


X71V2V'^ 

2 vl 


rrij 

7^ 


In 


In 


- 1 




+ 


o LL 


g^vf 


+ 9lvl{2 + 3ln^ 


- ^m%g^ ( - - In ^ 


(2 + 3 In + 2AmlY,^ (l - In - 2417?^^ In ^ 
+3mlY^ (^1 - In - SvlY* In + 8 [F{h, 1)]^ In ^ 




+4:mlG{h, 1) ^In ^ -1^+8 [F{H, 1)]^ In + 4mjjG{H, 1) ^In - 1^ 

1 ^A^ 2 fml^Vl ^ 3A71172 , A61171 A6117l^72^^ ^ \ 

+ -21^3-; 


6471^ 

+8 


m 


2R 

■12 


V 

1 - 1 ) “ n 2 A 5 i — 


2i7i 

A7in^ 


2173 


A7i17| 


Aei'^i Aei'ui'u^ 


^"V2i;i V viJ 4171 2171 2 ^ 

,0 2 f mj^Vi 1 3A71172 A61171 A 6117 i 17 ^A a 

+ 8 ">h* (— - - — + — - 


4172 

^H± 


mh 

In^ 


+4 


m 


2R 

12 


17i 

2 2 


17 \ 

1 - 1 ) “ n2A45 — 

Vv) } 


A7117I 


2i72 
A 71 I 72 


— Aeifi + 


In 

Aeiwii;^ 


17i 2 i7i 

^^wal { I -In^^ +5-1^2 f2 + 31n^) - 3m%g^ (7 - In ^ 




V 




2172 

1 


- 1 


In 


g? 


+ 


^ (2 + 3 in Y) + 24™^? (1 - l-Y) - in ^ 




Y4:rnf^G{h, 2) ^In —^ — 1^+8 [F[h, 2)] In —^ + Aw?^G(^H, 2) ^In 


+8 [F{H, 2 )fhi'^ 


m 


H 


- 1 
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Here 


G{h,l) 


G{h, 2 ) 


2 m?:? . 2 m??r;^ A7ir;2 3A6ifi\ 

5 o I 


. 2Ai- — + 

4 V r;ir;2 




V 2 J 


4Va 

X (^Ai + 


+2 ( Ai'Ui — 


( \ivl - \,_vi - Ivi - vl) + (Aei - An) _ ^ 

V viV 2 2 2 r;i 2 v 2 J 

'^12 (2 2'\ A7if2 SAsi^iA 

, r..2 ..2^ , , , , A7it'i , SAeir;?^ 


r’ir ;2 ^ 1^2 


2 iJ 

12 


+2 (A345r;2 + SAeifi)^ + OAei ( —m\2 + A345z;ir;2 + 


3A6ifi 


2\ 1 


+ 


Sa^G 


(hvi- X,vl - („5 - vi) + 3'‘’it-2(Ae.-A„) _ ^ ^ ^A 

V V1V2 2 2 z;i 2 r ;2 J 


X Ait>i - 


m 


2R 

12 


V 2 


2R 


3f2 (Aei — A 71 ) A 7 ii ;2 3A6iff 


+ 


4i;? 


+ 


4i;2 


+2 [ —m?? + A345 Wii;2 + 


3A6ifi 3A7 ii;| 


X (A345^^2 + 3A6lfl) 


1 2 m?? 2 mj^v‘^ 3 A 7 ii ;2 Aeii^i ^ 

4 V ^ V^V2 VnVl Vl V 2 J 


V 1 V 2 


+ 


V 2 V 1 


4\/a L 


A 2 X 2 "^1? / 2 2A , 3 i;ii ;2 .. , . A7ii;f AeiiifA 


/ m?? m?? 2 2 2 A 3 A 7 ii ;2 Aeii^iA 

(r ^ (”■ “ ”4 - J 

f X , " 1 ?? m?? 22 2 A , 3i;i (Aei - A 71 ) 3 A 7 ii;| Aei^iV 

^ +-J- 


+2 - 


+ 6 A 71 ( —m?? + A345i;ii;2 + 
1 


^ + ?^)+2(A2«..i+3A„„2)'' 


+ 


Sa^G 


X ' —X2V2 + 


m 


2 R 


fi 


12 , 7 2 21 "^12 , 3i;i (Aei - A 71 ) 3 A 7 ii;? Aeii;?^ 


+2 —4“ A345'Ui'U2 + 


3A6if? 3A7 ii;? 


1 2 


X (A345i;i + 3A7 ii;2) 


(A.5) 


(aix,; - A 2.1 - ^ - f|) + 

V V1V2 2 2i;i 2 i ;2 / 


(A. 6 ) 
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a{H,l) = 1 (2A1 - ^ + 

4 \ V1V2 V^V2 V2 ) 


+ 


1 


4^ 


(X,vt - X 2 vi -^^{vi- 4 ) + ^ 

V V1V2 2 2vi 2 v 2 J 


X Ai + 


m 


2R 

12 


V 1 V 2 


m 

v 


+2 ( Aiiii — 


m 


2R 

12 


V 2 


'-12 /,,2 ,,2A ^71 '^’2 I SAemA 
\v 2 ^ ^ 2 i;f 2 i ;2 / 
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